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I.  INTRODUCTION 


A.  Objecti ve 

The  purpose  of  this  report  is  to  describe  the  results  of  the  numerical 
computation  of  transonic  flow  over  an  ogive-cylinder-boattail  projectile  with 
a  simulated  rotating  band.  These  computations  are  compared  with  the  experi¬ 
mental  data  of  Miller.^  The  experimental  results  obtained  with  the  rotating 
band  are  also  summarized. 

B.  Background 


The  rotating  band  is  added  to  an  artillery  shell  in  order  to  impart  stabi¬ 
lizing  spin  to  the  shell  during  launch.  After  launching,  the  surface  protu¬ 
berance  of  the  band  adds  unwanted  drag  and  can,  through  changes  to  the  flow 
field,  cause  other  effects  on  the  aerodynamic  characteristics.  The  transonic 
region  is  particularly  important  because  most  aerodynamic  variables  (drag, 
pitching  moment,  side  force  and  Magnus  moment)  are  usually  at  their  maximum  in 
this  speed  range  and  thus  they  are  most  sensitive  to  changes  in  configuration. 
A  review  has  been  made  in  Reference  2  of  published  information  on  rotating 
band  aerodynamics  and  of  some  related  data  on  protuberances  and  steps.  It  is 
concluded  there  that  very  little  is  understood  about  the  magnitude  of  the 
effects  caused  by  the  addition  of  a  rotating  band  in  any  speed  range.  For 
example.  Figure  1  (which  is  taken  from  Reference  2)  shows  the  considerable 
uncertainty  regarding  the  rotating  band  drag.  Although  the  band  drag  may  only 
be  of  the  order  of  5  percent  of  the  total  drag  of  the  projectile,  the  data 
cited  in  that  figure  show  as  much  as  a  factor  of  two  uncertainty  in  the  rota¬ 
ting  band  drag  coefficient  at  transonic  speeds. 

The  primary  effect  of  flow  disturbances  caused  by  the  rotating  band  is  to 
increase  the  pressure  ahead  of  the  band.  This  is  because  of  the  turning  of  the 
viscous  layer  which  is  followed  by  a  series  of  expansion  regions  as  the  flow 
turns  to  move  over  the  band.  The  flow  then  recompresses  as  it  returns  to  the 
basic  shell  surface.  This  process  is  shown  schematically  in  Figure  2. 
Actually  some  of  these  turning  processes  are  caused  by  local  separation  and 
recirculation  regions  in  front  of  and  behind  the  band.  These  effects  have 
been  simulated  in  this  work  by  introduction  of  ramp  like  surfaces  ahead  of  and 
aft  of  the  band  in  order  to  produce  the  strong  pressure  disturbances  of  the 
protuberance.  If  this  can  be  shown  to  adequately  reproduce  the  effects  of  the 
band,  it  would  greatly  simplfy  the  prediction  of  the  aerodynamics  of  practical 
configurations. 

The  numerical  solution  technique  employed  here  is  the  azimuthal-invariant 
unsteady,  thin  layer  Navier-Stokes  code  which  has  been  developed  for 


3.  M.  C.  Miller,  'Vind  Tunnel  Measurements  of  the  Magnus  Induced  Surface 

Pressures  on  a  Spinning  Projectile  in  Transonic  Speed  Range,"  AIAA  Paper 
No.  83-1838,  July  1983. 

2.  J.  E.  Danberg,  'Numerical  Modeling  of  Rotating  Band  Flou  Field  and  Compar¬ 
ison  with  Eaperiment, "  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  Maryland,  ARBRL-TR-02505,  July  1983.  (AD  A13126Q) 
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application  to  projectile  shapes  by  Nietubicz,  Pulliam  and  Steger.^  The 
details  of  the  equations  used  and  of  the  method  of  obtaining  solutions  are 
described  in  the  following  sections  along  with  the  turbulence  model  employed. 

II.  EXPERIMENTAL  DATA 
A.  Facilities  and  Test  Conditions 


The  numerical  results  reported  here  are  compared  to  the  data  obtained  by 
Miller.^  The  main  emphasis  of  his  work  was  the  measurement  of  pressures  on  a 
spinning  projectile  configuration.  However, the  primary  interest  here  is  the 
nonspinning  case  so  that  comparison  with  the  numerical  calculation  is  possi¬ 
ble.  The  tests  were  performed  in  the  NASA  Ames  Research  Center  14-Foot  Tran¬ 
sonic  Wind  Tunnel  on  the  projectile  configuration  defined  in  Figure  3.  The 
conditions  of  these  tests  were  Mach  number  of  0.94,  a  Reynolds  number  of 
4  X  10^  1/Ft.  and  angle  of  attack  of  0,  4,  10°.  The  ogi ve-cyl inder-boattai  1 
model  was  instrumented  with  20  pressure  taps;  however,  only  three  of  these 
were  on  the  rotating  band.  Figure  4  shows  the  location  of  all  the  taps. 
Tests  were  performed  with  and  without  the  rotating  band. 

B.  Configuration  of  the  Rotating  Band 

The  rotating  band  model  afterbody  which  was  constructed  for  these  tests  is 
shown  in  Figure  5.  The  dimensions  were  based  on  measurements  obtained  from  a 
M549,  155  mm  shell  which  was  recovered  after  being  fired.  In  this  way  the 
band  was  made  to  simulate  the  conditions  in  flight  as  accurately  as  possible. 
The  crossection  of  the  band  is  shown  in  Figure  5.  Of  the  three  pressure  taps 
two  were  located  on  the  land  and  one  in  the  groove  of  the  band. 

III.  AZIMUTHAL-INVARIANT  THIN-LAYER  NAVIER-STOKES  CODE 
A.  Navi er-$tokes  Equations 

The  form  of  the  Navier-Stokes  solution  technique  employed  here  is  the  code 
developed  by  Nietubicz,  et  al^  and  applied  to  a  number  of  projectile 
configurations  as  reported  in  References  4  and  5.  The  axisymmetric  equations 


3.  C.  J.  NietubiaZj,  C.  J.  Pulliam  and  J.  L.  Steger^  "Numerica.1  Solution  of 
the  Azimuthal-Invariant  Thin  Layer  Navier-Stokes  Equations ^  "  U.  S.  Army 
Ballistic  Research  Laboratory ^  Aberdeen  Proving  Ground ^  Maryland,  AEBRL- 
TR-02227,  March  1980.  (AD  A085716) 

4.  C.  J.  Nietubicz,  "Navien-Stokes  Computations  for  Conventional  and  Hollou 
Projectile  Shapes  at  Transonic  Velocities,"  U.S.  Army  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-hlR-03184,  July  1982. 
(AD  A116866) 

5.  C.  J.  Nietubicz,  G.  R.  Inger  and  J.  E.  Danberg,  "A  Theoretical  and  Ea:peri- 
mental  Investigation  of  a  Transonic  Projectile  Flew  Field,"  U.S.  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland, 
ARBRL-MR-03291,  July  1983.  (AD  A131938) 
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can  be  written  in  transformed  body  oriented  coordinates  and  in  nondimensional 
variables  as: 


3q  +  3[.E  +  3G  +  H  =  Re  3  S 

where 


p 

p  u 

p  V 

p  w 


0 

0 

? 

A 

1 

+ 

II 

<  LU 

X 

S 

^ 

’  J  J 

X 

^z 
w  -  ^ 

T 

T 

The  countravariant  velocities  are: 
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and  the  ri-invariant  source  term  can  be  written: 
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where  R  is  the  radius  of  the  body  and  is  a  constant  of  proportionality 
between  the  circumferential  angle  and  the  transformed  coordinate  direction. 
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0.5(^2  +  72  +  w2)^  +  (a2)^[(Y  -  l)Prg] 
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These  equations  represent  the  time  average  turbulent  flow  equations  in  the 
following  sense: 

1.  The  velocities  denoted  with  a  tilde  are  mass-weighted  time  average 
quantities,  i.e. , 

u  =  'pu/p  (5) 


where  the  bar  superscript  refers  to  conventional  time  average. 
2.  The  effective  total  Prandtl  number  is: 


k  +  k  y/Pr  +  Ui./Pr 

1  (6) 

Cp(u  +  u  ) 


where  k^  and  y^  are  the  eddy  conductivity  and  eddy  viscosity,  respec¬ 
tively.  Inherent  in  this  form  of  the  equations  is  the  assumption  of  a 
Bousinesq  formula  for  the  time  average  Reynolds  stresses. 
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3.  The  thin  layer  approximation  has  been  applied  after  the  transformation 
to  body  oriented  coordinates  so  that  only  derivatives  with  respect  to 
the  "near"  normal  to  the  body  coordinate  (c)  have  been  retained  in  the 
viscous  terms. 

4,  The  consistent  equations  of  state  which  connect  "e  in  the  energy 
equation  with  the  mean  pressure  and  temperature  are: 


P  =  (y  -  1)[  e  -  0.5p(u2  +  +  w^)] 


T  =  p/{pR) 


(7) 


a^  =  yRT 


where  the  turbulent  kinetic  energy  has  been  neglected  relative  to  the 
mean  flow  kinetic  energy. 

In  the  present  non-spinning  application  of  these  equations,  the  body  sur¬ 
face  boundary  condition  becomes  the  adiabatic,  no-slip  condition.  Problems 
caused  by  the  singularity  on  the  axis  are  avoided  because  the  flux  vectors 
there  are  zero  by  symmetry.  Extrapolated  outflow  conditions  are  applied  at 
the  far  stream  boundary  which  includes  an  extended  sting  in  order  to  avoid 
computation  of  the  separated  near  wake  flow  field. 

The  equations  are  actually  solved  as  finite  difference  equations  written 
in  delta  form.  An  implicit  approximate  factorization  scheme  is  employed.^ 
The  algorithm  is  second  order  in  space  with  integration  of  the  unsteady 
equations  proceeding  until  steady  state  is  achieved. 


B.  Turbulence  Model 


A  two-layer  algebraic  turbulence  model  was 

these  calculations.  The  model  is  due  to  Baldwin 
boundary  layer,  the  model  can  be  described  by  the 


employed  to  determine  in 

and  Lomax. ^  For  an  attached 
following  equations: 


1.  Inner  layer 

o 

=  p£2  ja)| 


(8) 


6.  R.  Beam  and  R.  F.  Warming^  "An  Implicit  Factored  Scheme  for  the  Compressi¬ 
ble  Navier-Stokes  Equations,"  AIM  Paper  77-645,  June  1977. 

7,  B.  S,  Baldwin  and  H.  Lomax,  "Thin  Layer  Approximation  and  Algebraic  Model 
for  Separated  Flews,"  AIAA  Paper  78-257,  January  1978. 
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where 


l  =  Kz[l  -  exp(-z^/A^)] 


(z  =  normal  to  the  surface  in  the  untransformed  coordinate,  = 
u^z/v^.  A'*’  =  26).  |a)|  is  the  magnitude  of  the  vorticity  vector. 

2.  Outer  layer 


'^t  ~  ^cp  ^  ^wake  *^kleb 


(9) 


where  or  ^^/^max  smallest.  The 

quantities  and  Fj^gx  determined  from  the  function: 


F(z)  =  z|a)|  [1  -  exp(-z'^/A'^)]  (10) 

where  is  the  maximum  value  of  F(z)  and  z^^g^^  is  the  value  of  z  at 

which  it  occurs.  The  function  F|(ig|rj(z)  is  the  Klebanoff  interrnittency 
factor  given  by: 

''klebi")  =  -  5-5  r‘  (11) 

max 


The  quantity  U  is  the  maximum  velocity  in  the  profile  for  attached 
flow.  The  constants  not  already  defined  are:  K  =  0.0168,  C^p  =  1.6, 
Ckieb  “  0.3,  =  0.25  and  k  =  0.4. 

C.  Computational  Model  of  the  Rotating  Band 

Although  in  principle  the  Navier-Stokes  code  in  the  form  employed  here  is 
capable  of  computing  thin  regions  of  separated  flow,  the  method  is  far  from 
routine.  Thus  a  model  rotating  band  configuration  was  employed  which  provides 
attached  flow  over  a  protuberance  to  simulate  the  flow  above  the  dividing 
streamline;  however,  it  does  not  account  for  the  slip  boundary  condition  on 
that  streamline.  Figure  6  shows  the  mathematical  form  used  to  define  the  geom¬ 
etry.  The  height  of  the  band  is  h/D  =  0.0131.  The  length  of  the  quadratic 
ramps  ahead  of  and  behind  the  band  are  L  =  lOh  which  corresponds,  in  at  least 
a  rough  way,  with  the  extent  of  the  separated  region.  This  model  for  the  band 
is  identical  to  that  used  at  supersonic  speeds  as  discussed  in  Reference  2. 
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D.  Grid  Generation 


A  modular,  general  purpose  grid  generation  routine  as  described  in  Refer¬ 
ence  8  was  used  to  obtain  a  smoothly  varying  computational  mesh.  This  method 
allowed  arbitrary  body  geometry  and  investigator  selected  clustering  of  points 
near  the  band  and  other  critical  points  on  the  body.  Initial  computations 
were  performed  using  90  longitudinal  points  on  the  model  and  40  radial  points. 
The  number  of  points  available  on  the  rotating  band  was  very  limited  which 
reduced  the  resolution  in  the  region  of  rapidly  varying  flow.  A  new  grid  was 
selected,  employing  125  points  in  the  longitudinal  direction  with  100  on  the 
5.8  caliber  model.  Figure  7  shows  an  overall  view  of  the  grid  system  and 
Figure  8  shows  a  detail  near  the  rotating  band.  There  are  18  points  on  the 
band  with  4  points  each  on  the  forward  and  rear  ramps. 

For  the  ogive  nose  region  an  elliptic  grid  generator  has  been  employed. 
Straight  lines  normal  to  the  axis  of  symmetry  were  used  on  the  cylinder, 
rotating  band  and  afterbody.  An  exponential  stretching  of  points  away  from 
the  model  was  used  with  the  minimum  spacing  at  the  surface  being  0.00002  diam¬ 
eter.  The  computational  domain  extends  not  less  than  four  body  lengths  in 
all  directions. 

The  one-half  caliber  boattail  was  terminated  at  5.8  calibers  corresponding 
to  the  physical  model  length  and  a  constant  diameter  sting  was  continued  from 
that  point  to  the  downstream  boundary  of  the  computational  domain. 


IV.  RESULT 

The  numerical  solution  for  the  flow  over  the  projectile  permits  the  inves¬ 
tigation  of  many  details  of  the  flow  field  which  help  explain  the  relatively 
few  results  which  are  obtained  experimentally.  The  experiment  of  Miller  pro¬ 
vides  surface  pressure  data  at  20  points,  only  five  of  which  are  in  the 
vicinity  of  the  rotating  band.  On  the  other  hand,  the  numerical  solution  is 
limited  to  a  finite  special  resolution  which  affects  the  validity  of  the 
results.  Thus  the  comparison  of  experimental  and  computational  pressure  dis¬ 
tributions  are  used  to  establish  the  level  of  confidence  in  the  computation. 
The  examination  of  the  various  aspects  of  the  numerical  solution,  e.g.  Mach 
number  contours,  skin  friction  coefficients  and  velocity  profiles  gives  under¬ 
standing  of  the  physical  processes  involved. 

A.  Pressure  Coefficient  Distribution 


Figures  9  and  10  show  the  comparison  of  the  numerical  and  experimental 
surface  pressure  coefficient  distributions  obtained  at  Mach  number  0.94  with 
and  without  the  rotating  band.  Figure  9  shows  the  comparison  on  the  projec¬ 
tile  without  the  band  which  indicates  the  overall  capability  of  the  numerical 
technique  to  predict  the  pressures.  Figure  10  shows  the  comparison  with  the 


8.  J.  L.  Stegev,  C.  J.  Nietubias  and  K.  R.  Heavey,  "A  General  Curvilinear 
Program  for  Progeatile  Configurations,"  U.S,  Arrry  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-MR-03142,  October 
1981.  (AD  A1 07334) 
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rotating  band  on  the  shell  between  x/D  of  4.762  and  5.236  calibers.  There  is 
the  expected  rise  in  pressure  ahead  of  the  band  followed  by  the  expansion  as 
the  top  of  the  band  is  reached.  The  data  from  the  three  taps  on  top  of  the 
band  (Cp  =  -0.25)  are  in  good  agreement  with  the  calculation.  Behind  the 

band,  on  the  ramp  of  the  numerical  model,  the  computation  predicts  an 
extremely  sharp  expansion  with  a  minimum  Cp  of  about  -0.65,  followed  by  a 

recompression  to  a  pressure  corresponding  to  that  observed  ahead  of  the  boat- 
tail  without  the  band.  The  two  pressure  taps  located  in  the  region  are 
slightly  higher  than  that  obtained  numerically  but  of  generally  the  correct 
level.  From  the  start  of  the  boattail  at  5.3  calibers,  the  pressure  distribu¬ 
tion  is  somewhat  more  negative  with  the  band  than  that  obtained  without  it 
(see  Figure  9).  The  quality  of  agreement  between  theory  and  experiment  is 
essentially  the  same  for  both  configurations.  The  total  model  length  is  5.8 
calibers  and  the  relatively  poor  agreement  of  the  last  few  points  may  be 
attributed  to  the  sting  like  numerical  model  which  is  not  a  good  representa¬ 
tion  of  the  base  region  flow  field. 

The  comparison  between  computation  and  experiment  in  the  case  without  the 
rotating  band  (Figure  9)  shows  generally  good  agreement  with  some  discrepancy 
in  the  recompression  region  downstream  of  the  ogive-cylinder  junction  and  on 
the  boattail.  These  discrepancies,  of  the  order  of  8  percent  in  pressure  are 
probably  due  to  grid  resolution  of  the  expansion  and  subsequent  shock  wave. 
An  overall  improvement  in  the  representation  of  these  regions  was  obtained  by 
increasing  the  number  of  longitudinal  grid  points  from  90  to  125.  The  region 
near  the  base  is  poorly  predicted  again  because  of  the  sting  model. 

As  a  general  conclusion  it  may  be  said  that  the  experimental  pressure  dis¬ 
tribution  is  fairly  accurately  described  by  the  numerical  model  with  some 
qualifications  in  regions  of  minimum  pressure.  Unfortunately  no  pressure  taps 
were  located  in  the  expansion  region  just  behind  the  rotating  band  so  that  no 
conclusion  regarding  the  adequacy  of  the  flow  field  modeling  in  that  location 
is  possible. 

B.  Mach  Contours 


Figure  11  shows  the  calculated  Mach  contours  for  the  rotating  band  case. 
The  supersonic  pockets  just  after  the  ogive-cylinder  junction  and  behind  the 
rotating  band-start  of  the  boattail  are  clearly  evident.  Both  of  these 
regions  are  terminated  by  shock  waves  of  comparible  strength.  Apparently,  the 
subsonic  compression  region  ahead  of  the  band  involves  relatively  small 
increases  in  Mach  number  and  pressure  as  compared  to  the  effect  of  the  same 
surface  change  in  supersonic  flow.^’^  The  recompression  just  downstream  of 
the  rotating  band,  apparent  from  the  pressure  distribution,  does  not  show  up 
as  a  shock  wave  or  as  a  pocket  of  significantly  decreased  Mach  number.  The 
Mach  number  appears  to  remain  supersonic  throughout  the  region  until  the  boat- 
tail  recompression  shock  wave  just  ahead  of  the  base. 

C.  Skin  Friction  Coefficient 


Figure  12  shows  the  skin  friction  distribution  obtained  from  the  numerical 
solution.  Five  major  spikes  are  observed  related  to  five  regions  of  high 
pressure  expansion. 


14 


1.  Near  x/D  =  0.0;  this  is  the  expansion  created  by  the  hemispherical 
nose  cap-ogive  junction  which  forms  the  numerical  approximation  to  the 
blunt  nosed  projectile. 

2.  x/D  =  3;  the  junction  of  the  ogive  and  cylinder. 

3.  x/D  =  4.9;  at  the  leading  edge  of  the  top  of  the  rotating  band.  The 
local  minimum  ahead  of  this  peak  is  attributable  to  the  compression 
ramp  of  the  model  band. 

4.  x/D  =  5.1;  the  rear  of  the  band  which  is  followed  by  a  minimum  due  to 
the  recompression  on  the  ramp  and  the  short  downstream  cylindrical 
element. 

5.  x/D  =  5.3;  the  expansion  at  the  boattail  junction. 

One  significant  unexpected  result  is  the  negative  skin  friction  obtained 
on  the  boattail  which  suggests  separated  flow  at  the  boattail  shock  wave.  The 
configuration  without  the  rotating  band  shows  a  very  low  in  this  region  but 

no  negative  values.  This  shows  the  potentially  important  downstream  effect  of 
the  rotating  band  which  could  cause  some  change  in  the  boattail  drag  and  could 
explain  why  the  rotating  band  minimum  pressure  in  this  region  is  lower  than 
the  corresponding  value  for  the  no  rotating  band  case.  However,  since  the 
experimental  pressure  distribution  does  not  agree  with  the  prediction,  it  is 
uncertain  whether  or  not  separation  is  induced  by  the  band  in  the  experiment. 

D.  Velocity  Profiles 

Although  the  numerical  model  of  the  rotating  band  is  not  physically  cor¬ 
rect  because  it  does  not  properly  describe  the  separated  flow  regions  near  the 
surface,  the  strong  pressure  gradients  provide  a  severe  test  of  the  turbulence 
model.  Figures  13  and  14  show  the  development  of  the  velocity  profiles.  The 
curve  labeled  x/D  =  4.61  on  the  first  figure  corresponds  to  a  point  just  at 
the  begining  of  the  ramp.  Curves  labeled  4.76  and  4.85  are  on  the  compression 
surface  and  show  relatively  little  change.  Curve  4.92  is  on  the  top  of  the 
band  and  5.00  is  on  the  centerline.  Both  show  considerable  acceleration  of 
the  flow  in  the  viscous  layer  and  a  pronounced  overshoot  at  the  outer  edge  of 
the  boundary  layer.  Curves  5.04  and  5.08  of  Figure  14  continue  the  develop¬ 
ment  on  top  of  the  band  and  5.15,  5.24  and  5.22  are  on  the  downstream  ramp. 
Curve  5.15  is  just  off  the  top  of  the  band  and  shows  the  very  strong  local 
acceleration  experienced  by  the  flow.  The  subsequent  points  show  the  spread¬ 
ing  of  the  viscous  region  as  the  flow  begins  to  relax  back  toward  the  constant 
pressure  region  on  the  short  cylinder  before  the  boattail.  Consistent  with 
the  Mach  number  contour  plots,  the  regions  away  from  the  surface  show  an 
almost  monotonic  increase  in  velocity  as  the  flow  moves  over  the  rotating 
band.  The  boundary  layer  becomes  thinner  in  the  expansion  regions  with  strong 
overshoots  in  velocity  near  the  edge  of  the  layer. 

E.  Effect  of  Rotating  Band  on  Overall  Aerodynamic  Characteristics 

The  previous  sections  have  described  the  insight  gained  from  the  flow 
field  computation  at  zero  angle  of  attack  and  zero  spin  rate.  It  is  particu¬ 
larly  desirable  to  consider  the  effect  of  the  rotating  band  on  overall 
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aerodynamic  parameters  which  determine  the  shell  performance.  Although  the 

calculations  have  not  been  performed  at  angle  of  attack  or  with  spin  it  is 

possible  to  use  the  zero  angle  of  attack  calculations  to  help  explain  the 
experimental  results  obtained  from  the  NASA  Ames  tests.  Table  1  provides  a 
summary  of  the  integrated  data  obtained  with  and  without  the  rotating  band. 
Four  properties  are  considered:  normal  force,  pitching  moment,  side  force, 
and  Magnus  moment  all  in  terms  of  nondimensional  coefficients.  These  data 
were  obtained  by  integrating  the  measured  pressure  distributions  and  are 
broken  down  in  terms  of  the  contribution  from  the  ogive  nose,  the  cylindrical 
centerbody  and  the  afterbody.  The  rotating  band  pressure  distribution  is 
included  in  the  cylindrical  contribution.  The  normal  force  and  pitching 

moment  are  considered  to  show  the  effects  of  angle  of  attack  and  that  spin  has 

only  a  higher  order  effect  on  these  characteristics.  On  the  other  hand,  the 
side  force  and  moment  are  essentially  asymmetric  effects  created  by  the  spin 
of  the  projectile  at  angle  of  attack.  Figure  15  defines  the  coordinate  system 
employed  with  these  forces  and  moments. 

1*  Normal  Force  and  Moment.  The  computed  pressure  distribution  at  zero 
angle  of  attack  shows  that  the  pressure  increases  ahead  of  the  band  and 
decreases  behind  in  roughly  the  same  amount  so  that  relatively  little  direct 
effect  is  anticipated  on  normal  force.  However,  the  negative  pressure  coeffi¬ 
cient  region  is  considerably  reduced  and  thus  the  asymmetries  at  angle  of 
attack  which  give  rise  to  the  negative  boattail  lift  are  also  reduced.  This 
is  born  out  by  the  data  in  Table  1.  The  decrease  in  negative  lift  on  the 

boattail  in  the  amount  of  20  percent  produces  an  overall  increase  in  the  C.. 

N 

a 

of  about  3.6  percent  for  this  configuration.  The  moment  coefficient  is 

affected  in  a  similar  amount  but  opposite  sign.  Thus  a  decrease  in  the  nega¬ 

tive  lift  of  the  boattail  results  in  a  decrease  in  the  de-stabi lizing  moment 
in  the  net  amount  of  4.6  percent. 

2.  Magnus  Force  and  Moment.  Magnus  effects  may  be  expected  to  be  sensi¬ 
tive  to  large  gradients  in  the  surface  pressure  distribution  as  well  as  the 
magnitude  of  the  pressures.  The  addition,  the  rotating  band  creates  very 
severe  pressure  gradients  with  the  result  that  the  side  force  on  the  cylinder 
is  significantly  increased.  However,  the  decrease  in  pressure  on  the  boattail 
produces  a  compensating  effect  on  the  afterbody  so  that  the  overall  change  in 

side  force  is  of  the  order  of  1.7  percent.  On  the  other  hand,  the  location  of 

these  forces  produces  a  significant  13.5  percent  decrease  in  the  overall 
Magnus  moment  coefficient. 

As  an  additional  check  on  these  experimental  results,  the  data  are 
compared  in  Table  2  with  the  data  reported  by  Platou  and  Nielsen  in  Reference 
9.  The  and  results  from  integrating  the  pressure  distribution  are  in 

excellent  agreement  with  those  obtained  from  direct  wind  tunnel  force  and 
moment  balance  tests. 


9.  A.  S.  Platou  and  G.  I.  T.  Nielsen,  "Some  Aerodynamic  Characteristics  of 
the  Artillery  Projectile  XM549,"  U.S.  Army  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  Maryland,  ARBRL-MR-2284,  April  1973.  (AV 
910093L) 
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3.  Rotating  Band  Drag  Coefficient.  Under  the  assumption  that  the  numeri¬ 
cal  prediction  of  the  pressure  distribution  is  at  least  approximately  correct 
despite  the  limitations  in  the  model,  it  is  possible  to  estimate  the  drag  of 
the  band  itself  by  assuming  that  the  maximum  pressure  acts  on  the  forward  face 
of  the  band  and  that  the  minimum  pressure  acts  on  the  rear  face.  Thus  the 
change  in  the  drag  coefficient  due  to  the  band  is: 


)h/D 


which  in  this  case  is  approximately  ACj^  =  0.017.  Integration  of  the  pressure 
over  the  boattail  results  in  an  increase  in  Cq  in  the  amount: 

ACp,  =  0.005. 

^BT 


Giving  a  net  increase  in  drag  coefficient  of  0.022  or  about  9  percent  of  the 
total  drag  (assuming  Cq  =  0.25)  neglecting  any  change  in  the  base  drag. 


V.  SUMMARY  AND  CONCLUSIONS 

This  report  describes  the  results  of  applying  an  azimuthal-invariant 
Navier-Stokes  computational  technique  to  predict  the  flow  field  about  an 
artillery  shell  flying  at  Mach  number  of  0.94.  The  results  are  compared  to 
experimentally  measured  pressure  distributions  obtained  by  Miller^  in  the  NASA 
Ames ■ Research  Center  transonic  wind  tunnel. 

The  numerical  technique  is  essentially  that  developed  by  Nietubicz,  et  al^ 
which  includes  the  Baldwin-Lomax  turbulence  model.  In  order  to  avoid  comput¬ 
ing  separated  flow  regions  ahead  of  and  behind  the  rotating  band,  the  recircu¬ 
lating  regions  were  simulated  by  quadratic  ramps  of  length  equal  to  10  band 
heights  before  and  behind  the  actual  bands. 

The  presence  of  the  band  produces  a  significant  change  in  the  pressure 
distribution  with  a  high  pressure  region  on  the  upstream  ramp.  This  is  fol¬ 
lowed  by  an  extremely  short  expansion  zone  where  the  pressure  falls  to  a  level 
corresponding  to  supersonic  flow  on  top  of  the  band.  The  resulting  pressure 
agrees  with  those  obtained  experimentally  although  the  computation  did  not 
include  the  details  of  the  lands  and  grooves  of  the  experimental  rotating  band 
configuration. 

The  pressure  tap  locations  in  the  experiment  were  not  placed  so  as  to  con¬ 
firm  or  contradict  the  existence  of  a  second  sharp  expansion  which  occurs  just 
downstream  of  the  rotating  band.  The  numerical  prediction  shows  a  drop  in 
pressure  coefficient  immediately  behind  the  band,  followed  by  a  recompression 
on  the  aft  ramp  to  about  the  same  level  as  on  top  of  the  band  and  consistent 
with  the  pressure  obtained  just  ahead  of  the  cyl inder-boattail  junction  with¬ 
out  the  band.  The  pressure  on  the  boattail  appears  to  be  only  slightly 
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modified,  but  small  changes  can  have  significant  effects  on  the  projectile 
drag.  The  boattail  contributes  an  important  fraction  of  the  total  drag  at 
transonic  speeds. 

The  reason  for  the  decreased  pressure  on  the  boattail  is  attributed  to 
flow  separation  at  the  boattail  shock  wave.  This  separation  region  is 
observed  in  the  skin  friction  distribution  which  is  not  observed  in  the  calcu¬ 
lations  without  the  band. 

The  experimental  results  regarding  the  overall  aerodynamic  coefficients  of 
the  shell  can  be  interpreted  by  considering  the  numerical  solution  even  though 
the  current  numerical  results  apply  only  to  zero  angle  of  attack  and  zero 
spin.  The  change  in  boattail  pressure  distribution  because  of  the  rotating 
band  induced  separation,  results  in  slightly  increased  normal  force  coeffi¬ 
cients  at  angle  of  attack  and  a  decrease  in  pitching  moment.  This  is  consist¬ 
ent  with  a  reduced  negative  lift  on  the  boattail  attributed  to  separation. 

The  band  induced  change  in  the  side  force  is  negligible  but  this  is  due  to 
a  compensating  decrease  in  boattail  contribution  matching  an  increase  on  the 
cylinder  because  of  the  rotating  band.  The  Magnus  moment  is  the  most  signifi¬ 
cantly  affected  quantity.  The  experimental  data  indicate  a  13.5  percent 
decrease  because  of  the  band. 

The  numerical  pressure  data  were  used  to  estimate  the  change  in  drag  coef¬ 
ficient  caused  by  the  rotating  band.  The  estimated  Cq  is  increased  by  slight¬ 
ly  less  than  7  percent  due  to  the  pressure  acting  directly  on  the  band.  The 
decreased  pressure  on  the  boattail  also  increases  the  drag  by  about  2  percent 
so  that  the  total  increase  in  drag  is  9  percent  under  the  conditions  consider¬ 
ed  in  this  report. 
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Figure  1.  Rotating  Band  Drag  Coefficient 
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Figure  2.  Rotating  Band  Flow  Field 
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Figure  3.  Conical  Boattail  Model  Configuration  Showing  Dimensions 
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Figure  4.  Pressure  Tap  Locations 
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Figure  5.  Rotating  Band  Configuration 
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Rotating  Band  Numerical  Model 
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Figure  8. 


Detail  of  Grid  in  Vicinity  of  Rotating  Band 
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Figure  9,  Pressure  Coefficient  Distribution  on  Model  without 
Rotating  Band  and  Comparison  with  Experiment 
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Figure  10.  Pressure  Coefficient  Distribution  on  Model  with  Rotating  Band  and  Comparison  with  Experiment 
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Figure  11.  Mach  Number  Contours 
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Figure  12.  Skin  Friction  Coefficient  Distribution,  Rotating  Band  Model 
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Figure  13.  Velocity  Profiles  in  the  Vicinity  of  the  Rotating  Band  (Forward  Half  of  Band) 
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Figure  14.  Velocity  Profiles  in  Vicinity  of  the  Rotating  Band  (Aft  Half  of  Band) 
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Table  1.  Aerodynamic  Coefficients  Based  on  Integration 
of  the  Measured  Pressure  Distribution 


MODEL  CONFIGURATION: 


8  INCH  DIAMETER  MODEL 
3  CALIBER  OGIVE 
2  CALIBER  CYLINDER 
5  CAL  BOATTAIL 
REF  eg,  AT  Z/L  -  .62B 


TEST  CONDITION: 

MACH  .94 
a  =  10° 

R^  =  4x10®/FT 
p  =  -4900  RPM 

p  =  -.162 


TERM 

C|y  (OGIVE) 

a 

Cm  (CYLINDER) 
a 

Cm  (BOATTAIL) 

a 

Cm  (total) 

a 

C,^  (OGIVE) 

C  (CYLINDER) 
"(1 

C^  (BOATTAIL) 

C,„  (TOTAL) 

o 

^cp/L 


Cy  (OGIVE) 

P 

Cy  (CYLINDER) 
P 

Cy  (BOATTAIL) 
P 

Cy  (TOTAL) 

P 

C,^  (OGIVE) 

P 

C  (CYLINDER) 

V 

C  (BOATTAIL) 

^  P 

C^^  (TOTAL) 
^cp/L  (MAGNUS) 


ROTATING  BAND 
OFF 

1.94 

.60 

-.31 

2.24 

3.47 

-.36 

.60 

3.71 

.33 

-.004 

-.092 

-.080 

-.176 

-.003 

.063 

.148 

.208 

.836 


ROTATING  BAND 
ON 

1,95 

.62 

-.25 

2.32 

3.47 

-.39 

.46 

3.54 

.35 

-.004 

-.149 

-.026 

-.179 

-.003 

-.139 

.043 

.180 

.804 
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Table  2 


Comparison  of  Aerodynamic  Coefficients  from  Two  Experiments 


TERM 

Cy^  (TOTAL) 
Cn  (TOTAL) 

Cy^  (TOTAL) 

C,  (TOTAL) 

np 


INTEGRATION  OF  DIRECT  FORCE  AND 


SURFACE  PRESSURE 

DATA 

MOMENT  DATA.  1 
BRLMR  2284 

a  =  10° 

-.179 

-.175 

.180 

.180 

P 

II 

0 

-.104 

-.090 

.090 

.085 

Test  Conditions 


M„  =  0.94 

GO 

pd/2V  =  0.164 
Rotating  Band  On 
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LIST  OF  SYMBOLS 


time  average  of  the  velocity  of  sound  squared 

drag  coefficient  (based  on  projectile  cross  sectional  area) 

pressure  coefficient 

specific  heat  at  constant  pressure 

derivative  of  the  normal  force  coefficient  with  respect  to  angle  of 
attack 

derivative  of  Magnus  moment  coefficient  with  respect  to  spin  rate 

derivative  of  pitching  moment  coefficient  with  respect  to  angle  of 
attack 

derivative  of  side  force  coefficient  with  respect  to  spin  rate 
body  diameter 

rotating  band  diameter  in  calibers 
energy 


E,G,H,S  flux  vectors  of  transformed  Navier-Stokes  equation 
h  height  of  rotating  band 

J  transformation  Jacobian 

k  molecular  thermal  conductivity 

k|.  turbulent  eddy  thermal  conductivity 

L  body  length 

i  Prandtl  mixing  length 

M^  free  stream  Mach  number 

Pr  Prandtl  number 

Prg  effective  Prandtl  number  (see  Equation  6) 

Pr^  turbulent  Prandtl  number  =  l<^/(CpP^) 

p  time  average  pressure 
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LIST  OF  SYMBOLS  (Cont'd) 


q  flux  vector  of  dependent  variables  in  the  transformed  Navier-Stokes 

equati ons 

R  radius  of  the  body  in  n-invariant  Navier-Stokes  equations 

^  gas  constant 

r  radius  of  the  body 

r|3  radius  of  body  with  rotating  band 

U,V,W  contravari ant  velocity  components 

u,v,w  mass-weighted  time  average  velocity  components 

u  shear  velocity 

Greek  Symbols 
a  angle  of  attack 

Y  ratio  of  specific  heats 

K  universal  constant  in  Prandtl  mixing  length  equation 

p  molecular  viscosity 

turbulent  eddy  viscosity 

5,n,c  transformation  coordinates  in  axial,  circumferential  and  radial 
di recti ons 

p  time  averaged  density 

T  transformed  time 

(j)  circumferential  angle 

kinematic  viscosity  evaluated  at  the  body  surface 
ci)  vorticity 

Subscri pts 


RB 

rotating  band 

BT 

boattai 1 

t 

turbulent 
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LIST  OF  SYMBOLS  (Cont'd) 


e  effective 

"  free  stream  conditions 

Superscripts 
_  time  average 

~  mass-weighted  time  average  (see  Equation  5) 
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